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Abstract 

We consider a diffusion model with limit cycle reaction functions. In an unbounded domain, 
diffusion spreads pattern outwards from the source. Convection adds instability to the reaction- 
diffusion system. The result of this instability is a readiness to create pattern. We choose the 
Lambda-Omega reaction functions for their simple limit cycle. We carry out a transformation 
of the dependent variables into polar form. From this we consider the initiation of pattern to 
approximate a travelling wave. We carry out numerical experiments to test our analysis. These 
confirm the premise of the analysis, that the initiation can be modelled by a travelling wave. 
Furthermore, the analysis produces a good estimate of the numerical results. Most significantly, 
we confirm that the pattern consists of two different types. 
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I. INTRODUCTION 



Morphological patterning, such as animal coat markings, may be caused by a chemical 
field {ll,^,!^. The Turing model now has strong experimental support [4, for a review], jsl,^]. 
This mechanism has some limitations, and so we continue to investigate variations of the 
model [7]. 

n 

We consider the standard reaction-diffusion model with the addition of convection [8|, 
y, [l^. The more general form is a system with advection, as in fll^. Our motivation is 
theoretical, to see the effect of convection on the robustness of the pattern formation. We 
consider that a weak amount of convection, or a similar effect, may be present in vivo. 
This may soften the standard requirements for the formation of pattern, increasing the 
applicability of the model. 

There may be a direct biological application, such as the formation of the vertebral 
precursors (somitogenesis) . The organism growth could produce such a convective effect 12| . 
However, convection-induced patterning is not generally supposed for somitogenesis fl3|. 

Our system does not fit the standard Turing analysis since our reaction functions are 
already unstable. However, we have shown previously that an unstable function may produce 
a Turing pattern [14]. Furthermore, we have shown that, in the presence of convection, we 
can see similar behaviour irrespective of the stability of the fixed point 14 1. 

Experimentally, convection is introduced from the boundaries 16]. The effect of 



boundary conditions is also likely to be relevant in vivo. The boundary can have a significant 
effect on the pattern formed [9]. However, we follow the suggestion of Cross and Hohenberg 
to examine the system in a boundary- free environment first [17] . 

In numerical simulations we make form an initial point disturbance. Pattern is formed, 
spreading outwards from its initiation point. The disturbance is oscillatory and complex. 
Our focus here is the way in which the pattern propagates. 

The speed of propagation has been given for a two-species system with no convection pi 
and one with equal convection [9]. A related theoretical study considers the relationship 
between the onset of the instability and the longer-term behaviour of the non-equilibrium 
state [isl. Here the simplest reaction-diffusion system with oscillatory kinetics was consid- 
ered, and complex behaviours were found. 

Limit cycles are inherent in these oscillatory models, since an unstable spiral at the 
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steady state must be bounded [20[]. In this paper, we seek to exploit the hmit cycle to aid 
our analysis. To this end, we select the Lambda-Omega reaction functions. These are chosen 
to produce the simplest limit cycle possible: the unit circle. 

The reaction functions chosen in a full model correspond to actual reaction mechanisms. 
Examples of these include the Belousov-Zhabotinsky reaction, and Brusselator-type systems 
such as the chlorine dioxide-iodine-malonic acid reaction. The Schnakenberg model can be 
considered an intermediate step towards these functions. For an initial investigation, we have 
chosen simple, related functions. The intention is to extend the analysis to more realistic 
models. Secondarily, the results in this simple case may form a useful basis for comparison 
to more complex systems. 

We transform the problem so that the oscillations are removed, or at least reduced. This 
brings the onset of the pattern into much sharper relief. The two aspects - the pattern, 
and the initiation of the pattern - are clearly distinguished. Thus this onset of pattern is 
effectively a phase transition. 

The transformation is to convert the dependent variables into polar form. We then 
consider that the formation of the limit cycle approximates a travelling wave, and so we 
employ a Fisher solution to the problem [2V\. The relative sizes of the parameters affects 
our analysis: we consider various situations. We carry out numerical experiments to discover 
which of these estimates are valid. 

II. LIMIT CYCLES AND THE X-oj FUNCTION 

We suppose that some form of chemical mechanism is the underlying basis for biological 
pattern formation. If the chemicals are well-mixed, then the law of mass action valid, and 
an ODE is an appropriate model: 



We consider reaction functions chosen for their simplicity in the form of the limit cycle 
produced: 




(1) 



g = u + v{l — V? — f ^) . 



(2) 
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FIG. 1: Phase space for the hmit-cycle reaction ([I]), ([2D, given by numerical solution. The phase 
curves spiral out from the steady state (the origin) to meet the limit cycle (the unit circle). The 
trajectories starting far away from the steady state spiral into the limit cycle. Solid curves starting 
with circles are the trajectories. 

This is a simple form of the X-uu class of functions. Here the steady state is the origin 
and the limit cycle is the unit circle. There is no parameter to determine the stability: the 
steady state is unstable and the limit cycle stable. This behaviour is clear in the {u^ v) phase 
plane (see Figure [1]). The circled points on the diagram show the start of the phase plane 
trajectories. The phase curves spiral out from the steady state to join the limit cycle. The 
trajectories starting far away from the steady state spiral into the limit cycle. 

These functions are chemically unrealistic, as they stand. However, we can relate these 
functions to the Schnakenberg reactions, chosen to be the simplest chemical form which 
can produce a limit cycle. This mechanism was proposed theoretically, but has been used 
as a model for actual reaction mechanisms 
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23[ |. The key reaction there is the cubic 
autocatalytic one, U + 2V 3V. Using the law of mass action, this produces the term uv'^. 
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The reverse step gives v^. These types of terms form the core of the X-uu functions, the cubic 
terms. 

We consider a new coordinate system for the dependent variables ix, ^' as a polar form r, 9 
as follows: 

tan^ = - . (3) 
u 

By differentiating these identities, we transform the first order ODE ([T]) into: 

r r(l — r^) , 

e = 1, (4) 

The 9 equation clearly decouples and is resolvable. The remaining single-species ODE in r 
has steady states at —1,0, 1. The negative state is unrealistic because r must be positive, 
zero is unstable and one is stable. The point r = 1 is our limit cycle: the unit circle. 
Furthermore, 9 ^ time and so progression around the limit cycle is constant. 

Thi s ty pe of analysis is often used in similar cases, such as the analysis of travelling wave 



trains [2J]. 



III. PATTERN FORMATION 



Biological pattern formation is by definition spatially differentiated. There are many 
models for this. One is that chemicals diffuse within an organism, then cells respond differ- 
ently dependent on the concentration of one of these chemicals. If the concentration of the 
chemicals has formed a pattern, then this is refiected by the cells. 

For an initial analysis, we choose the simplest case: one spatial dimension and some dif- 
fusion. By diffusion we refer to the averaged gross effect of random motion of the chemicals; 
with passive movement this is equivalent to Brownian motion. In the case of general reaction 
functions, this system is the one proposed by Turing [Jb]. 

We add a convective term to the system. This can be considered a disturbance to the 
system, as an investigation of stability. For a small amount of convection, this could cor- 
respond to axial growth [12]. However, for convection to be significant in axial growth we 
would require extremely slow diffusion of the chemicals to fit the model. 
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The general system is then as foUows: 

ut = eiu^^ -pu^ + f , 

vt = £2Vii -qv^ + g ■ (5) 
We consider different convection on each species (p ^ q). This effect could be created in a 

n 

chemical flow reactor, with one of the reactants held flxed in a packed bed p^. For biological 
applications, one of the chemicals may be held within a cell, while the other flows freely. If 
both chemicals flow freely, it is possible that their movement may be hampered by obstacles 
such as the extra-cellular matrix. In this case, a larger molecule may be affected more 
strongly than a smaller one, and different convection speeds may result for each chemical. 
This concept parallels that of different diffusion rates, as in the Turing model. 

We consider that p > 0. If this is not the case, we make the transformation {p^q^i) 
(— p, — g, — ^), giving a positive value for p with no change to the form of ([5]). We can remove 
one of the convective terms by a simple change of coordinates x = ^ — pt: 

Ut = SiU^x + / , 

Vt = ^2Vxx ^ (6) 

with / and 5^ as in (Ej) and Si^ £2 and 7 positive constants. From the original system ([5]) we 
have ^ = q — p. If 7 turns out to be negative, we can reverse the sign as we did for p and 
q. Equivalent ly, we choose the coordinate change x = pt — ^ previously, and set 7 = p — g. 
In the case where the convection is the same on both species {p = q)^ then both convection 
terms are removed (7 = 0). This recovers the basic Turing model. 

We have reaction functions that have limit-cycle behaviour and convection, which is 
known to drive instability. The appearance of pattern is then to be expected, although 
the form might be more difficult to predict. Satnoianu, Merkin and Scott [25|] studied a 
similar system to ([6j) previously, with Schnakenberg reaction functions in place of the X-u 
ones. They found that periodic behaviour is emergent in the system over a broad parameter 
range. 

In the numerical experiment, we start at the steady state (0,0), except for a small dis- 
turbance at X = 0. We try to simulate a boundless environment - to this end we find zero 
derivative boundary conditions the most effective. The initial disturbance propagates and 
becomes pronounced, forming a regular pattern with aligned oscillations. The propagation 
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FIG. 2: Pattern found for a diffusion system with convection and limit-cycle reaction kinetics 
The initial disturbance propagates and becomes pronounced, forming a regular pattern with aligned 
oscillations. The propagation is linear, forming a V-shape. The convective effect is slight: the 
pattern is skewed slightly to the right. In this case, the direction of the internal oscillations is 
distinct from the angle of propagation of the pattern. This is a numerical solution of ([6]) using 
NAG D03PCF, plotting species u with 7 = = ^2 = 1- The reactants are initially at steady state: 
(u^v) = (0,0), with a small disturbance at x = 0. The boundaries are held at zero derivative: 
Ux = 0^ Vx = 0. 

is linear, forming a V-shape. The convective effect is to skew the pattern to the right (Fig- 
ure [2j). Given the parameters in the figure (7 = ^1 = ^2 = 1), the direction of the internal 
oscillations is distinct from the angle of propagation of the pattern. The emergence of this 
pattern is the main focus of our study. 
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IV. TRAVELLING WAVE ANALYSIS 



Before we start any specific analysis, we rescale the spatial variable x to remove one of 
the parameters. We choose x = s/siy. This yields the system: 

vt = Vyy + evyy -^Vy + g , (7) 

where 7 = ^/s/si. Following suite, we expect £2 = ^2/^1, but we go one step further in 
defining s = £2 — I. This produces a symmetric first spatial derivative on both equations, 
with £ the difference between the two diffusion rates. 

We wish to again convert u^v into the polar form (r, 0). During the conversion, we see 
that there is grouping of the terms parameterised by £ and 7: 

n = Tyy - + -{eVyy " y\ + T {\ " T^) , 

T U 

Ot = Oyy + 2j-dy + -[SV yy " ^ V y] + 1 . (S) 

We complete the transformation: 

n = (1 + £Sin^6') Tyy 

+ (2e sin 9 cos 9.9y — 7 sin^ 9) Vy 

+ (1 - - 6'^ + £ sin 6' cos 9.9yy - e sin^ 9.9l - 7 sin ^ cos 9.9y) r , 

9t = (1 +£COS^6') 9yy 

+ (2^ + 2/-^ cos^ ^ - £ sin ^ cos 9.9y - 7 cos^ 9^ 9y 

^l^^_IjnLZJl]L^i^ecos9 . (9) 
r 

If we translate our numerical results into this polar form, the onset of the pattern becomes 
very clearly demarked: the new coordinate r transitions sharply from to 1. The primary 
behaviour is that of a travelling wave, propagating outwards. The internal behaviour of the 
pattern, the steady oscillation, is reduced to a secondary effect. There are minor ripples in 
the established solution: the limit cycle settled on by the PDE is not the unit circle, but is 
close. There is a difference in behaviour between the left and right sides of the pattern: the 
angle of alignment and the frequency of the pattern is different on either side (FIGURE [3]). 
Next we examine this travelling wave analytically. 
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FIG. 3: Polar form of the pattern. The new coordinate r transitions sharply from to 1. This is 
a travelling wave, propagating outwards. There are minor ripples in the established solution: the 
limit cycle settled on by the PDE is not the unit circle, but is close. The oscillation is at an angle 
to the travelling wave front. There is a difference in behaviour between the left and right sides of 
the pattern: the angle of alignment and the frequency of the pattern is different on either side. 
There is also a clear centre to the propagation, roughly at x = 7/2t. This is a numerical solution 
of ([6j) using NAG D03PCF, plotting species u with = £2 = 1, 7 = 2. The reactants are initially 
at steady state: (ix, v) — (0, 0), with a small disturbance at x = 0. The boundaries are held at zero 
derivative: Ux — Vx — 0. 

A. Simple version 

We consider the very simplest situation: equal diffusion on both chemical species {e = 0) 
combined with no convection of v (7=0). We make a simplifying assumption for an initial 
analysis: 9y ^ 0^ 9yy ^ 0. The problem reduces to 



n Tyy + r(l - r^) 



(10) 
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Then 9 decouples, as for the ODE. The solution for ^ is ^ = t, to within an arbitrary 
constant. 

The remaining equation in r is of the form similar to the Fisher equation and should 
therefore yield a propagating wave solution. This analysis is covered in greater detail in I2J]. 
We look for a solution of the form 

r{y,t) = R{z), z = y-ct. (11) 

which gives 

R'' + cR' + R{l-R^) =0 . (12) 

We carry out a phase plane analysis in the (i?, phase plane. We find that c > gives a 
stable point at [R^R') = (0,0), which suggests the solution we are looking for. The point 
i? = 1, which corresponds to our limit cycle, is always a saddle point. We are looking for a 
phase plane trajectory that leaves this saddle point and goes to the zero steady state: this 
will be our travelling wave solution. 

Small wave speed, c < 2, gives a spiral in the phase plane, which would imply r < 
at some point on the trajectory. From the definition of r we know this to be impossible 
and thus this is unrealistic for a travelling wave solution. For c > 2 we have a node and 
the trajectory discussed above, leaving from the saddle point will head directly to the node, 
remaining in the fourth quadrant of the phase plane and therefore realistic. This trajectory 
equates to the travelling wave that we are looking for. We expect that the least speed wave 
will be achieved so we predict a wave of speed c = 2. Converting this back into our unsealed 
system, we have cr = 1-s/e{ to the right. There is also a solution to the left: = —2^Je{. 
In the original system, we have wavespeed = p± 2y^. 

B. Full system 

We now apply the method for the full system. Here we may have different diflFusion on 
the two species {s ^ 0) or some convection (7 > 0). We again make the assumption 9^ ^ 0, 
9^^ ^ 0. We introduce the travelling wave coordinate z = y — ct and look for a solution 
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R{z) = r{y,t) . We also linearise the system, dropping the term: 

(1 + esin^ 6) R" + (c - Tsin^ 6) R' + R = , 
eR" - ^R' 

^ ' sin^cos^ + 1 = 0. (13) 
R 

Using the same approach as before, we find c = 7 sin 6:sin^. Converting this 

back to the unsealed coordinate system we have c = 7 sin^ 0±2^ Si cos^ + £2 sin^ 9. We see 
that 7 and £2 remain hnked by the sin^ term. These parameters occur on the v differential 
equation, and v = rsinO. 

In the original system we have wavespeed = pcos^ 9 + qsin^ 9 ± 2\^£i cos^ 9 + £2 sin^ 9. 
Again, the parameters for each species u and v are linked by the functions cos^ and sin^ 
respectively. 

This general result for the wavespeeds gives a dependence on the angle variable 9. We 
consider some general numerical experiments to assess the behaviour of the system. What 
we find is two distinct types of behaviour. We refine our wavespeed estimates in light of 
these particular cases. 



C. Small parameter behaviour 

For parameters close to those of the simple system, namely small convection (7 <C 1) and 
near-equal diffusion (6:1 ^ 6:2), we expect the behaviour to remain similar to that of the simple 
system. Previously we have seen ODE-like and Turing-type behaviour in a convection-free 
system [14]. In both cases, we see that the internal angle of the pattern (the group direction) 
is not in line with the wavefront. As such, the oscillation is perceptible at the wavefront. 
This suggests that an average value approximation of 9 may be an appropriate estimate. 

In this case we approximate sin^ 9 with its average value: 1/2. This yields c ^ 7/2 ± 
2^1 + £/2. In the unsealed coordinate system we have c ^ ± \/2^j£\ + ^2- 

If we consider the limit as the amount of diffusion tends to zero, we see that the centreline 
of the propagation is at 7/2. In this situation, the convection does not increase the separation 
of the left and right wavefronts; the spread of the propagation is solely due to diffusion. In 

v^q. ^ ^ 

the original system the incident fronts are at wavespeed ^ ± v 2v^i+^- 

For dimensional reasons, it is appropriate to write £i = £\ . Then the diffusive contribu- 
tion to the wavespeed becomes \f2\/ + £2 ^ which is reminiscent of Pythagorean form. 
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More simply, the term is additive and the wavespeed is proportional to the diffusion of each 
species. 

D. Large parameter behaviour 

As the strength of convection 7 increases, we see an alignment of the incident wave with 
the internal group direction. This means that the value of 6 is approximately constant at 
the onset of the wave. We investigate the effect of this on the wavespeed estimate. 

We consider the range of values possible for the right wavespeed with a fixed value for 9. 
Since sin^ has a range of [0, 1] then c must range at least between 2 and 7 + 2\/l + e. This 
second value is not necessarily greater than the first, for example 7 = and ^ = — 1/2 gives 
\/2 1.4, which is less than 2. However, once 7 > 2 the second value is greater than the 
first, so we take this as an estimate for the transition to large parameter behaviour. Thus 
the range of possible values for the right wavespeed is cr G [2, 7 + 2^1 + 

On the left we have a similar result: Cl G [—2, 7 — 2^1 + ^] is the range once 7 > 2. We 
note that for sin^ = 0, ^ = (mod tt) and for sin^ = 1, = 7r/2 (mod tt). 

We propose that the maximum spread of the instability is achieved. That is, the 
wavespeeds are maximised within their given possible ranges. For this, the left wavespeed 
must reach its lower limit, and the right its highest. This g ives C]^ — 2 and Cr 
7 + 2^1 + Then the onset angle to the left is expected to be approximately (mod tt), 
and to the right of the right wave the angle should be roughly 7r/2 (mod tt). 

We see this is the case in Figure [H There is a smooth transition between the left and 
the right of the pattern. The oscillations at the centre line (approximately 7/2) are roughly 
perpendicular to the centre line itself. 

In the unsealed coordinate system, for strong convection, we have Cl ^ and 
Cr ^ 7+2-y/^. In the original system, this is left wavespeed ^ p— 2y^ and right wavespeed = 
q + for q > p. In the other case p > q^ there is a reversal and left wavespeed ^ q — 
and right wavespeed = p + 2^Je{. 

This analytical prediction was made from observing that the internal waves become 
aligned with the onset of the instability. However, when we make the assumption that 
the spread of the propagation is maximised, it follows that the onset angle Q become con- 
stant. Furthermore, the suggestion of 7 = 2 as a transition value between behaviours occurs 
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FIG. 4: Large parameter behaviour of 9. The oscihations are ahgned with the travehing waves. To 
the left of the left travelling wave (the light, solid line), the angle is roughly 0, and to the right of 
the right wave (the dark, dashed line) the angle is approximately 7r/2. There is a smooth transition 
between the left and the right of the pattern. The oscillations at the centre line (7/2; medium, 
dot-dashed line) are roughly perpendicular to the centre line itself. This is a numerical solution 
of (P) using NAG D03PCF, plotting species u with 7 = 5, £1 = 62 = 1. The reactants are initially 
at steady state: (ix, v) — (0, 0), with a small disturbance at x = 0. The boundaries are held at zero 
derivative: Ux — Vx — 0. 

naturally from the algebra. 

V. NUMERICAL EXPERIMENTS 

Having carried out an analysis together with some initial numerical investigation, we now 
conduct a more comprehensive experiment. We vary our parameters over a wide range and 
measure the incident wavespeeds in each case. The results confirm our analyses, as we can 
see in Figure [51 

For the left wavespeed the difference between the low and high parameter estimate is 
pronounced, so the transitional behaviour is noticable. There is less difference between the 
estimates on the right, and the transition is close to the intersection of the estimates. This 
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(a) Left wavespeed. The darker surface is the low (b) Right wavespeed. The darker surface is the 
parameter estimate: 7 — \ple\ the hghter surface low parameter estimate: 7 + \ple\ the lighter 
is the high parameter estimate: —2. Both surface is the high parameter estimate: 

estimates hold well with a smooth transition 7 + 2y^. Both estimates are good and there is 
between behaviours. a neat transition between behaviours. 

FIG. 5: Numerical data compared to analytical estimates. The estimates correspond to the data 
over the full range. The low parameter estimate holds for longer for high e than for low. The 
points are data from numerical solutions of ^ using NAG D03PCF, plotting species u. The 
parameter e — 82/ Si'^ in this experiment Si — 1 so e — 82- The reactants are initially at steady 
state: {u^ v) = (0, 0), with a small disturbance at x = 0. The boundaries are held at zero derivative: 
Ux = 0^ Vx = 0. 

seems quite remarkable, that the behavioural transition occurs v^here the estimates intersect. 
As a result, there is no discernable area of transitional behaviour. 

In both cases the lov^ parameter behaviour endures longer for s high, compared to s low. 
This suggests that the transition, caused primarily by the convection, 7, is held in check 
more strongly by the diffusion on v^here the convection is applied. 

VI. ABSOLUTE VERSUS CONVECTIVE INSTABILITY 

There is a concept of stability in systems such as these, v^hich considers v^hether the 
system at a particular spatial point, once destabilised, v^ill return to the steady state. If so, 
the system is classified as convectively unstable. If the chemical species do not return to 

14 



stasis then the system is deemed absolutely unstable. 

In line with our analysis, we define absolute instability as sign {left wavespeed) ^ 
sign {right wavespeed)^ with equivalence giving convective instability. Combining this with 
our results, we see that for small parameters we require 7^ > 8 (^1 + £2) for convective in- 
stability in the unsealed system. This inequality is unlikely to be satisfied whilst keeping 
within the small parameter regime. 

However, if we consider the original system, the condition is {p + g)^ > 8 (^1 + ^2), which 
is quite possible to satisfy while still remaining in the 'small parameter' regime). For exam- 
ple, p = q = 2 gives 7 = 0, and then Si = £2 = I satisfies the condition. 

For large parameters we always have absolute instability in the unsealed system. Referring 
to the observed data (Figure [5j), we see the left and right wavespeeds always have diflFerent 
signs. Thus, we can classify this system as only exhibiting absolutely unstable behaviour. 

However, if we consider the original system ([5]) we find that q > p > or p > q > 
will give convective instability rather than absolute instability. 

We see that, in the original system, we can always shift the stability from absolute to 
convective with sufficient convection applied to the system. In the ffow coordinate system, 
the stability is essentially absolute. However, when we convert to this system we make an 
arbitrary choice to remove the parameter p rather than, for example, q. Then the distinction 
between absolute and convective instability becomes similarly arbitrary. 

VII. AN ALTERNATIVE PERSPECTIVE 

At the start of our analysis, in section [TTIl we chose a new coordinate system for the 
purpose of eliminating one of the convective terms. This reduction in the number of param- 
eters proved useful in the ensuing work. However, the choice of coordinates, as we observed 
previously, was arbitrary. 

Now we consider an alternative coordinate system, one centred on the convection. 

From ([5]) we choose x = ^ — t. This gives the system as: 

vt = ^2^xx -jvx + g^ (14) 
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q-p 

where 7 = — ^ — . We can recalculate the apparent wavespeeds from this perspective. For 
our 'small parameter' regime we have c = ±\/2V^i + ^2- Here we see that we have aligned 
our system centrally with the flow. 

In our large parameter case, we flrst consider q > p. This ensures 7 > 0. Now the left 
wavespeed is = — (7 + 2y^) , and the right wavespeed is c^^ = 7 + 

p-q 

In the other case p < q^ we deflne 7 = > 0: 

vt = e2V^x+lVx + g . (15) 

Now the wavespeeds are = — (7 + 2^^), and the right wavespeed is = 7 + 

Consider the right wavespeed. In each case it is a function of only one of the diffusion 
parameters. We examine the equation for the relevant species. In both cases, we flnd a 
negative convection term for that species. This suggests that negative convection drives the 
onset of instability. 

For the left wavespeed the connection is with the positive convection. However, in this 
direction x is decreasing. From this perspective, the signs of the convective terms change. 
In this sense then, the connection is the same. 

If we understand that the negative convection drives the instability wavefront, then we 
can consider only the associated species. For example, looking for the right wavespeed 
in (fT5l) . we see that u is the relevant species. We make the assumption that the other 
species v is much smaller than and so disregard it. Then we can carry out a travelling 
wave analysis for this single species. The result is exactly the same as we found in our polar- 
form analysis, with a much simpler approach. We can conflrm this behaviour by observing 
that the initiation angles and 7r/2 correspond to u and v only, respectively. 

We can go one step further in unravelling the puzzle. If we simply separate the species 
- perhaps we could call this a super-linearisation - we can produce two wavespeeds, one 
from each species. In our example (fT5l) . for the right wavespeed we have the possibility of 
7 + 2y^ from the u equation and —7 + 2y^ from v. Then we apply the 'maximisation' 
principle suggested earlier. For large 7, we can be sure that the wavespeed given from the u 
equation will maximise the propagation. Thus we can flnd the wavespeeds by applying two 
simple concepts to the problem. 
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The behaviour for smaU convection rehes on some more complex way of combining the 
two species. The wavespeed in this case is dependent on the strength of diffusion of both 
the species. Therefore the polar- form approach given previously retains some merit. 



VIII. DISCUSSION 

The addition of convection to the reaction-diffusion model has produced many patterns j^, 
261]. Advection, being the broader definition, has potential to produce even more. Aspects 



of the pattern have been observed, such as the different behaviour to the right and the left 



of the free-forming pattern 



25|. 



Previously analyses have considered a variety of models, most restricted in different ways. 
No movement for one species [27], equal diffusion and equal convection [9] and fixed diffusion- 



convection ratio 



281] have all been considered. Here, we have considered the most general 



diffusion and convection possible, as in more recent work [26]. 

Much previous work was devoted to the discussion of 'flow and diffusion distributed 
structures' (29|. This is spatial pattern with no variation in time. Our reduction of the 
system to two parameters allows for a clear and concise analysis of the system. However, 
after our change to convective coordinates, the meaning of 'time-constant' is lost. In fact, 
in our reduced system, a time-constant solution is only possible for the trivial case, 7 = 0. 

The standard analysis to date has been a linear existence analysis for pattern, involving 
dispersion relations. This follows the Turing derivation, but is difficult for this more complex 
model. In our analysis we have employed a different approach. We have not directly consid- 
ered the internal behaviour of the pattern, but instead looked for approximate descriptions 
of the propagation speed of the pattern. This information has already been found for other 
versions of this system (less general diffusion or convection, more complex reactions) |^, Q. 

To understand the understand a system well, it is best to isolate influences [17]. This is 
why we consider the pattern forming away from the spatial boundaries. In this instance we 
also have symmetrical, parameterless reactions, giving the most elementary form possible. 
Previous work has shown boundary effects mixed in with other behaviour [261] , and that the 
boundary can control the existence of pattern [d]. 

Our analysis is the natural one given the system and its behaviour. We have combined two 
simple techniques to reduce a complex problem to a more manageable one. We consider only 
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the (apparently) linear aspect. Our analysis predicted the speed of the incident travelling 
waves over the full parameter range. In combination with the numerical investigation, we 
demonstrate two different behaviours of the system. There is a distinct transition between 
them, with the transition region indicated by the analysis. Our work also confirms the 
validity of the Fisher-type analysis in this two-species system. Furthermore, we give some 
insight to the different behaviour to right and left of the pattern. 

Having successfully analysed this simple X-uu system we must attempt a similar analysis 
on a more complex system. The next step is to consider a more realistic set of reaction 
functions, such as the Schnakenberg; ones. In this case, we already know that first-order 

n 

effects will come into play [30]. This should lead to a clearer understanding of the Turing 
phenomenon. 
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